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Integrability in general tetrad formalisms is reviewed, following and clarifying work 
of Papapetrou and Edgar. The integrability conditions are (combinations of) the 
Bianchi equations and their consequences. The introduction of additional con- 
straints is considered. Recent results on the conservation of constraints in the 1+3 
covariant formulation of cosmology are shown to follow from the Bianchi equations. 



1 Introduction 

In a number of recent works^'Bifl calculations have been made to show that 
when the governing equations of relativistic cosmology, in the covariant 1+3 
formalism, are divided into evolution equations and constraints, the constraints 
are preserved under the evolution. Unfortunately these attempts do not wholly 
agree with one another, though all agree the constraints are conserved. 

On general principle, the only integrability equations that should arise are 
the Bianchi identities (regarded as equations rather than identities, in tetrad 
formalisms) and their consequences, and the derivatives of any additional con- 
straints imposed (e.g. the perfect fluid form). Thus the cited results on the 
evolution of constraints should follow immediately from the Bianchi equations. 
Here we show explicitly that this is true, noting some subtleties involved. 

Integrability for general lelxad formalisms has been discussed by Edgafl, 
following work of PapapetroTB'Lia. The first section of this paper reviews and in 
part clarifies that treatment. Section 2 considers the introduction of additional 
constraints and section 3 shows how the conservation of the first two of the 
cosmological constraints previously discussed can be derived from the Bianchi 
equations, and reconciles this with previous treatments. A fuller discussion 
of these issues, in particular completing the discussion of the 1+3 constraints, 
will appear elsewhere. 
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2 Integrability in tetrad formalisms 

In the Cartan approach to curvature, there are two basic structure equations 
for a set of basis one- forms a;°. For a Riemannian space these are|^ 

du;" -r% A u;'' (1) 

whicli utilises the connection one-forms T"'b, and 

dT\ + T\AT'h = R\ (2) 

which defines the curvature 2- forms i2°h. The integrability conditions of these 
equations arc simply the conditions d^a;" = and d^F^f, — 0, the first and 
second Bianchi identities, the first of which can also be considered to be the 
Jacobi identities for a triple of the basis vectors dual to the w", 

J^fccd = ([eh, [ec, ej] + [e^, [e^, Cf,]] + [e<j, [e,,, 6^]]) = . (3) 

Written in terms of forms, these Bianchi identities read 

A , (4) 
dR^b - R^c A r"fc + r"^ AR'^b^O ■ (5) 

Note that (||) and (^) guarantee the integrability of (||). 

Tetrad formulations normally start (though they need not) with a basis in 
which the metric is fixed, i.e. if gab is the scalar product of and e^, 

dgab = d(ea.efc) = . (6) 

Tetrad indices can be lowered and raised with gab and its inverse. Together 
with which can also be written as the commutator equations 

[ea,eb] = -2r''[ab]ec, (7) 

and fixes the skew parts T^'^[bc] of the coefficients in F^b = T°-bci^'^, @ deter- 
mines r°f,c- Thus, under the condition (^), (|^) is equivalent to the definition 

Vab = Lu%eb'ea''-^ = -io%.,eb'ea'' (8) 

in terms of the cooordinate components (indices i, j, k, . . .) of the basis vectors 
and one-forms. Papapetrou and Edgar both remark on this equivalence, but 
without explicitly emphasizing that it only holds when (g) is assumed. 



"Index and sign conventions here follow Kramer et alu unless otherwise stated. 
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For a Ricmannian space, the components R'^bcd of the curvature two-form 



R\ = ^R\adOj''Aoj^ (9) 

must satisfy 

Rated ~ —Rbacd = — Rabdc ~ Rcdab, Ra[bcd] — . (10) 

In a tetrad formulation, it is usual to regard gab as given and the com- 
ponents Cfc*, V^ab and R°'bcd as functions to be solved for. The equations for 
them are then the tetrad component forms of (0), (H) and (|^). Since (H) is 
now regarded as an equation, one has to adjoin its integrability condition to 
close the system. Note that it is assumed that the components V^ab satisfy 
the symmetry condition Vcab + Tacb = and that the R°'bcd satisfy (|To|): both 
these are conditions which are implicit in the usual discussions, but they are 
important and should be made explicit. 

The equations of the tetrad formalism can thus (always assuming the given 
constant metric and the above symmetry conditions) be written as 

= 5T\, = T\ab] - ^^e^b'^a]. , (11) 

= SRabcd — Rabcd — ^abd.c + ^abcd — bd^ aec + ^"^bc^aed — '^^^ [cd]^ abe ,{^'2) 
= SBabcde = Rabcd-e + Rabde;c + Rabec-d ■ (13) 

Here a comma denotes a partial or directional derivative, and semi-colon de- 
notes the covariant derivative for which the V^bc are the connection coefficients 
(and which is the metric connection if (11) and (^) hold); this differs from 



Edgar's use of the semi-colon for a directional derivative of a tetrad compo- 
nent. The Bianchi identity ( p^ has of course to be expanded in terms of 
directional derivatives and the V^bc in order to give a set of equations entirely 
in terms of the declared set of variables. 



The quantity SRabcd in (g^ does not obey ( |10[ ) though it does satisfy 



SRabcd = —SRbacd = —SRabdc ■ (14) 

In 4 dimensions, there are six distinct skew pairs of indices so a tensor obeying 
( |r^ has 36 distinct components. For SRabcd these 36 can be decomposed 
into a Riemann tensor part with the symmetries (^) and 20 independent 
components which can be further split into the Weyl and Ricci tensor parts 
SCabcd and SRab, a part SRabcd - SRcdab (15 components) and a part SR^abcd] 
(1 component). The last two parts give SJ°'bcd- 

These equations are not exactly Cartan's, due to substitutions and the 
taking of tetrad components, so the integrability conditions are combinations of 
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the Bianchi equations. From Edgar's calculations, transcribed to my notation, 
these integrability conditions are 

= sj'^bcd - sv^iM.c] - <5r°[(,|e|(5r''crf] - sr^^iji.dr'^'^d] + st^^^iJ"'''^^ 

+^r"[6|e|r''cd] + SVelc^'^M] + [cd^'^ b]e , (15) 
= SBabcde — {5Rabcd)-e — {SRabde):c — {SRabec);d 

^6r„fc[c|^g|^r^(ie] - STab^SJgcde , (16) 
— r]^'^'^^ {{SBabcde) j — SRabgc^J^ def — R^ [a\cdRg\b]ef ~ SRabcd,gST^ ef -{^7) 

One can check that no further conditions arise. 

These results show that there is some redundancy in (pl])-(p^, e.g. that 
if (pl])-(p^ are satisfied so are (p^ )-(p^) and (p^. However, it is important, 
as Edgar has noted, that (p^)-(|l3|) are not a complete set of equations. This, 
and misunderstandings about the role of and (^), have contributed to un- 
certainty about whether the tetrad equations form a complete sjistem and if 
so which equations are required (segJor example^ Chandrasekhaii3) . The sit- 
uation has been discussed by EdgafiO. 

One could take a different set of variables. For instance, one could adjoin 
fimctions supposed to be the values of Rabcd;e', then (^, or more precisely the 
definition of Rabcd-e in terms of Rabcd and V^bc, are algebraic equations for 
Rabcd;e, and ( p7| ) is a differential equation, giving further integrability condi- 
tions. Finally (as Dr. A. Rendall has reminded me) one should note that in 
general integrability conditions are necessary but not sufficient to guarantee 
existence and uniqueness of solutions, which depend also on analytic properties 
of the equations and data. 

3 Additional constraints 

The equations above assume a completely general spacetime and choice of 
tetrad. Additional equations and integrability conditions arise if one imposes 
additional conditions which amount to replacing some of the inequalities which 
must hold for the general case by equalities. Among these, for equations of 
state which limit the energy-momentum tensor, are Einstein's field equations. 
Edgar reminds us that the above tetrad equations remain consistent in the 
case of vacuum in general relativity. For the cosmological case this is also true 
for (for instance) perfect fluids. 

One could distinguish between two different cases of choice of equality, or 
set of equalities, as follows. 

''This reference was drawn to my attention by Prof. B. Edgar. 
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First, one could choose a tetrad in a way that is compatible with the general 
case under study. The derivatives of this constraint will then give additional 
conditions. However, the processes of imposing the constraint and solving the 
equations must commute; that is to say one could in principle have solved the 
equations in a general basis and then chosen the tetrad by applying a Lorentz 
transformation at each point (i.e. algebraically). Therefore the extra equations 
in this case must be just the differential equations necessarily satisfied by a 
tetrad chosen after solving the original tetrad equations. In the cosmological 
context an example of a condition of this type is given (assuming the energy- 
momentum obeys the usual energy conditions pwhich guarantee the existence 
of a timelike eigenvector of the Rifjci. tensor) by the requirement of zero heat 
flux, Qa — in the usual notatioitatJ. 

A second type of equality or set of equalities that can be imposed is one 
which gives a consistent specialization of the general system of equations but 
which is not generally possible just by choice of tetrad: for example, the con- 
ditions that a spacetime is vacuum, of Petrov type D, and has shearfree and 
geodesic principal null congruences. These can, at least intuitively, be regarded 
as specifications of subspaces, in some suitable function space, such that solu- 
tions of the tetrad equations which start in the subspaces lie wholly in them. 
In cosmology, the assumption that the matter content is irrotational 'dust' 
(with no acceleration) is such a set of conditions. 

One often wants to work out such a set starting from some incomplete set 
of conditions. An example is given by the investigations, in cosmology, aimed 
at deciding what solutions, if any, satisfy conditi^ris such as the vanishing of 
the electric or magnetic parts of the Weyl tensooEJtj. This can become very 
confusing when the conditions are substituted into the general tetrad equations 
and integrability conditions are then considered. It seems to me, though I have 
not so far found a proof, that one can always arrive at the consequences in a 
clearer way by taking the extra conditions, say = for some set of indices A, 
and just adding the equations C'^^ab... = until the system closes, rather than 
intermixing the new conditions with the original equations before constructing 
the integrability conditions for them. 



4 Application to the 1+3 covariant formalism in cosmology 

The approach referred to here is the on based on a timelike congru- 
ence witlLunit tangent vector u": this is a 'threading' rather than 'slicing' 
approaclo. It is covariant in the sense that it is coordinate-independent, and 
invariant if is invariantly-defined. The metric is taken to have signature +2 
and one defines hab = gab + UaUb- The kinematic quantities 8 (expansion), aab 
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(shear), ojab (vorticity) and m° (acceleration) are then defined by 

Ua-b = ^Qhab + <Jab + UJab - UaUb (18) 

where 

CTab = (Tba, Cr^a = 0, (Tafcu'' = 0, LUab = ^tOba, UJabU^ = . (19) 

From the point of view above, one can regard this as an incomplete or- 
thonormal tetrad formalism in which the spatial triad has not been speci- 
fied (but in principle could be, if w° is invariantly defined, by following the 
methodology of the approach to characterizatiaa. of spacetimes due to Car- 
tan and Karlhede: see e.g. MacCallum and Skeai3). In this way of looking at 
the method, the kinematic quantities are (some of) the connection coefficients 

i be- 

The energy-momentum tensor is split by 

Gab = Rab - \Rgab = Tab = tJi-UaUb + ^q^^a^b) + PKb + T^ab (20) 

where qaU°- = ,Tiab = T^ba-, a = and TTahU^ = 0. The Weyl tensor Cabcd 
can be written in terms of electric and magnetic parts, Eab and i/afc, which are 
both symmetric traceless tensors orthogonal to u", as 

Cabcd — cjacEbd — QadEbc — QbcEad + QbdEac 

+r]abeUcHd'' - llabeUdH/ + TJcdeUaHb" - TjcdeUbHa" (21) 

where gab = gab + "^UaUb and rjabc = Vabcdu'^ is a three-dimensional volume 
element (in the tangent space at a point in the general case where ujab 7^ and 
hence the planes defined by hab are not surface-forming). Later I will refer to 
(combinations of) the equations (p^), etc., as S cia a nd so on. 

One now proceeds to split the equations (|l2| - |l3|) , as far as possible, by 
invariant projections parallel and perpendicular to u'^, ignoring those equations 
which cannot be written in terms of the above variables. The equations are 
separated into evolution equations, giving values of derivatives Q — Q-au"" of 
quantities Q, and constraints which do not involve such derivatives. Note that 
for a tensor projected orthogonal to u° the projection of Q differs from Q only 
by multiples of and ii^ with Q, so when studying conservation it is sufficient 
to study the projection. 

In the general case, the evolution equations consist of three Jacobi identi- 
ties giving the evolution of uJab = ^abc^'^, a combination of the Einstein equa- 
tions, the Raychaudhuri equation, giving 8, five components of SRabcd giving 
(Tah, and 14 Bianchi equations giving /i, qa, Eab and Hab- The constraints con- 
sist of one Jacobi identity giving uj°'-^a, three Einstein equations giving qa, five 
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components of SRabcd giving Hat and 6 components of SBabcde giving D°'Eab 
and D°'Hab, where the operator Z?" is defined for any tensor by 

DaT'- ''e...f = ha'h^ . . . h\h/ . . . h /TP-" r...s-b • (22) 

The general version of the tetrad equations is not needed in full in this pa- 
per, and so is omitted. Note that the evolution of some quantities remains 
undetermined and depends on, for example, an equation of state of matter. 

We note that in this system, compared with a general tetrad system, we 
have omitted 12 Jacobi equations and 6 components of SRabcd- The latter are 
those field equations, in a tetrad whose spacelike basis vectors are indicated 
by use of Greek sub- and superscripts, which require -Rq,/3. Without these, of 
course, the system is incomplete, and in particular the aab equations require a 
combination ol Eab and -Kab rather than these two, as might seem more natural, 
being separated. When one can make a separation, i.e. when one has a complete 
tetrad, it is usual to distinguish the field equations Gap = Ta/3 which give (Tq^ 
and equations for Eab and GabU°'u^ which are constraints. This organization 
of the possible linear combinations of the equations is especially appropriate 
when considering a Cauchy problem, since the constraints GabU^ = TabU^ are 
conserved by the evolution equations Gap = Tap as a consequence of the 
contracted Bianchi identities, as is well-known. 

The incompleteness of the_system of equations without Gap — Tap or an 
equivalent is discussed by Elli£3 and VeldeiQ. In the latter work, however, the 
missing equations themselves are not written down (they are only implicit in 
the statement in words of the meaning of Rap) and the linear combinations 
taken give dab in terms of Eab (only, among the components of Rabcd) and 
constraints for Rap- This system is then treated as a complete system of 
evolution and constraint equations, although it is still missing the equations 
giving Rap in terms of the connection or metric and is perhaps therefore better 
regarded, in the context of the study of evolution (though not in some of 
Velden's other discussions) as just adjoining to the usual system, with the 
equation SRap — omitted, some equations giving algebraic definitions of 
additional quantities which would, if Rap were in fact calculated, be that 
tensor. In a complete tetrad formalism, the 12 extra Jacobi identities include 
the evolution equations for the remaining V^bc insofar as these are fixed without 
conditions on the tetrad. 

It is now my aim to show that the discussion of evolution of constraints 
in the various works cited can be more simply derived, and more deeply un- 
derstood, by considering specialization of the conditions dl^p^. Here I will 
illustrate this only for the irrotational dust case, and only for two of the con- 
straints discussed in previous works; a fuller account will appear elsewhere. 



7 



For definiteness, I shall follow Maartenfl, though the other related papers use 
similar methods and notation. From the present point of view, the method is 
to assume that SV^bc = = SJatcd = SiTab = Sp and that SRabcd = = SBabcde 
for the components giving evolution equations. In the equations one can set 
u'^ = = u'^ = q'^ — TTab = P, which in particular means all Jacobi identities 
are ignored. For reasons which will become apparent, however, I will retain a 
non-zero for the moment. The evolution equations are then 

- SRoo = e + + + aaba^'' , 

= S{Eab - ^TTab) = ha'^hb &cd + fQcTab + 0'c<aO'b>'' + Eab , 

~ (5Bian(/i) = /i + 6^ , 

- (5Bian(<7a) = ha^'q' + aabq' + ^Qqa , 

= 6Bian{Eab) = ha'^h'^Ecd + QEab - icrc<aEb>'' - V(a\cd\D''Hb)'^ ~ ^H(Tab , 

= 6Bian{Hab) = ha^h'^Hcd + QH^b - ia,<aHb>^ + V{a\cd\D^E^'' , 

where I have adopted the convention that angle brackets denote the traceless 
symmetric part. The names 6Roo and so on have been introduced to make 
the following discussion easier to read without referring back to numbered 
equations. 

The constraints are 

= Sqa = qa + D^dab - \DaQ , 
= 5Hab = Hab - VialcdlD'^'Tb)^ , 

= ^Bian(i^^K6) = D'Eab - fi^aM " Vabc^' dH^" - ^aabq' + ^Qqa , 
= 5Bian{D''Hab) = D^'Hab + Vabc^rKiE"' - ^Vcd(a^"b)q'' ■ 

From (5Bian((7a) the evolution equation for dqa should be 

= ha''6q, + CTabSq'' + leSqa . (23) 

However, to caJculate, from the above system, the evolution of Sqa (denoted 
by Maartenfl), one has to eliminate 8 and &ab using SRqq — and 5{Eab — 
^TTab) = 0, which, from the Bianchi identities, naturally brings in (5Bian(D''i?afc) 
(the precise combination is given by a component of Thence we should 

have 

= ha'^hb + ^abh'' + i^Sqa + 5Bmn{D'' Eab) , (24) 

Since we are assuming the V^bc are correct, the initially zero qa has to 
be given the value —6qa while calculating the evolution of 6qa, and similarly 
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the initial Hab has to be replaced by Hat ~ SHab {SHab is denoted — C^^ by 
Maartenfl). Thus if = D'^Eab - \DaH - Vabc'y^ d.H'"^ , then 5Bieav{D^Eab) = 
Cl + ■nabcC^dSH'"^ + \aab5q'' - I'^ha leading to 

= hJ'Sqb + fflahV + QSqa +Cl- VabcCT^C^^'' , (25) 

This is precisely Velden's result. Maartenfl gives 

= ha'Sq, + eSqa +Cl- 2rjabcT\C^''' , (26) 

The discrepancy arises because the modified Ricci identities with non-zero 
SRabcd have only been used in some of the places where they should appear. 
Since the V^bc are taken to be correct, the SRabcd terms have to be added to 
the commutators wherever the corresponding Rabcd appear. For example, one 
should use 

(D'Sab)- = D'Sab - lOD'Sab ^ <j'"D,Sab + Vabc{H\ + C''\)S^'' - §5^^^ , 

n (27) 

rather than the similar equation with Sq^'' = = C^'^ used by Maartenfl. 
Applying the same method to 6 Hab, ^Bian(iJa6) would give 

= ha'hb^'SH.d + eSHab - 3a,<aSHb>' ~ hcdi,aO\)bq'' , (28) 

(in which the possible bEbc term is taken as zero due to i'Kab — and S{Eab — 
iTT^b) = 0). Thus 

= Cl, + eel, - ia.^aCl^" + \iicdia<y\)C^'' ■ (29) 

Maartenfl has 

= Cl, + eel, + vcdiaCT^jC'" . (30) 

This discrepancy is again due to an inconsistent way of taking a commutator, 
which should in this context read, in Maartens' notation in which curl Sab — 

VcdiaD'^Sb)", 

(curl S'ab)'^ curl Safe - iScurlS'ab - Cre^VcdiaD" Sb)'^ 

+ 3{Hc<a + Cl<a)^b>'' + ^Vcd(aS''b)C^'^ , (31) 

Velden has 

- Cl, + eel, + hcdi.ao^p" , (32) 

but the source of this discrepancy is not known to me. 

It is intended to give the remaining calculations for the irrotational dust 
case, and the results for more general cases, elsewhere. 
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